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ABSTRACT 

We present resolved-star VI photometry of the Local Group dwarf irregular galaxy IC1613 reaching 
I ~ 23.5, obtained with the wide-field camera at the 2.5 m Isaac Newton Telescope. A fit to the stellar 
density distribution shows an exponential profile of scale length 2.'9 ±0.1 and gives a central surface 
brightness jivfi = 22.7 ± 0.6. The significant number of red giant branch (RGB) stars present in the 
outer part of our images (r > 16.'5) indicates that the galaxy is actually more extended than previously 
estimated. A comparison of the color-magnitude diagrams (CMDs) as a function of galactocentric 
distance shows a clear gradient in the age of its population, the scale length increasing with age, while 
we find no evidence of a metallicity gradient from the width of the RGB. We present quantitative 
results of the recent star formation history from a synthetic CMD analysis using lAC-STAR. We find 
a mean star formation rate of (1.6 ± 0.8) x 10"'^ Mq yr~^ kpc~^ in the central r < 2'.5 for the last 
300 Myr. 

Subject headings: galaxies: dwarf — galaxies: individual (IC1613) — galaxies: irregular — galaxies: 
stellar content — galaxies: structure — Local Group 



1. INTRODUCTION 

IC1613 is a typical dwarf irregular galaxy concerning 
its luminosity, metallicity and star formation rate (SFR). 
In fact, it serves as the prototype for the DDO type Ir V. 
It is a low surface-bright ness galaxy with a moderate 
luminosity {My = —15.3. Ivan den B;erghll2000 ) located 
730±20 kpc from our Sun. This distance c orresponds to a 
distance modulus (m-M)o = 24.31±0.06 polphin et all 
I2001h and a scale of 3.54. PC per arcs econd. The recent 
calculation oflPietrzvhs ki et al.l (|2(20fD from near-IR pho- 
tometry of cepheids gives (m - M)o = 24.291 ± 0.035, 
for which the authors claim a total uncertainty of less 
than 3%. While the spatial extent estimate of lAblesI 
gives an optical size of 16' x 20', the obs ervation 
of car bon stars up to 15' from the center by Alber t et al.l 
(j2000) indicate that it is virtually twice that size. Its 
high galactic latitude, in the southern hemisphere, con- 
fers it low extinction and color excess. Her e we adopt 
a reddening of E{B -V) ^ 0.02 ± 0.02 from lCole etall 
([l999). 

The H II reg i on rn etalli city of IC1613 was mea- 
sured by iTalentI (|1980l see ISkillman et"aLl 119891) to be 

12 4-log(0/H) = 7.86 from spectrop hotometric obser- 
vations of the [OIII]A4363 line, while iLee erall ()2003f ) 
found 12 + log(0/H) ~ 7.62. This corresponds to 
[Fe/H] = -0.8 and [Fe/H] = -1.07 dex respectively 
(jSkillman et"al][2003 . hereinafter S03). The mean color 
(F-/)_3.5 of the RGB at Mj = -3.5 gives [Fe /H] =-1.3 
for th e old and intermediate-age population l|Lee et al.l 
Il993f ). This low overall metallicity and the high gas con- 
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tent (iHoffman et al.|[l996l ) suggest a primitive state in its 
evolution. 

The star formatio n his tory (SFH) has been studied 
quite extensively by [SOl for an HST/WFPC2 field lo- 
cated 7.'4 kpc southwest from the center. Their conclu- 
sion is a relatively constant SFR over a long period, with 
the oldest population being more than 10 Gyr old. They 
also summarize the results about structure and stellar 
content of the whole gal a xy fro m the literature. More 
recently, iBorissova et al.l ()2004D analyzed 60 OB asso- 
ciations, apparently correlated with H I I regi ons stud- 
ied kinema tically bv Lozinska va et al.l (|2003l, see also 
ISilich et al.ll2006f) . fMagrini et all (|2005D announced the 
detection of two candidate planetary nebulae. 

In this paper we describe a wide- field survey of IC 16 13: 
an overview of the observations and data reduction is pre- 
sented in section 121 and the resulting color-magnitude di- 
agram (CMD) is described in iJ31 In 21 we examine the 
extent and morphology of the galaxy, as well as the spa- 
tial structure of the different populations in the CMDs. 
An analysis of the recent SFH at different galactocen- 
tric radii is given in ^ Finally, in fJHwe summarize the 
results and present our conclusions. 



2. OBSERVATIONS & DATA REDUCTION 

Observations of IC1613 in Harris V and Sloan Gunn i' 
were conducted on 5 nights between November 1999 and 
September 2000 using the Wide Field Camera (WFC) at 
the 2.5 m Isaac Newton Telescope (INT) of the Observa- 
_ torio del Roque de los Muchachos. The WFC is a mosaic 
cpptgj|£.||j^ made up of four 2048 x 4096 CCDs, with a pixel 
size of 0.33". The total field of view is about 34' x 34', 
covering most of the galaxy. The total integration times 
in V and i' were 3660 and 1830 seconds, respectively. A 
detailed observing log is presented in Table [TJ 

Overscan trimming, bias subtraction and flat-field cor- 
rections were performed using the standard routines in 
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TABLE 1 
Journal of Observations. 



UTC Date 


Time (UT) 


Filter 


Exp. Time (s) 


Air Mass 


1999 Nov 06 


22:23 


V 


60 


1.14 


1999 Nov 06 


22:32 


V 


1200 


1.13 


1999 Nov 06 


22:53 


V 


1200 


1.12 


1999 Nov 06 


23:14 


V 


1200 


1.12 


2000 Aug 10 


04:38 


I 


30 


1.12 


2000 Aug 10 


04:41 


I 


600 


1.12 


2000 Aug 10 


04:54 


I 


600 


1.12 


2000 Aug 10 


05:05 


I 


600 


1.12 



10 arcmin 



IRAF'^ . The i' images have also been corrected for fring- 
ing effects. The DAOPHOT II/ALLS TAR and ALL- 
FRAME programs ()Stetsonlll98^ . Il994f ) were then used 
to obtain the instrumental photometry of the resolved 
stars from the four individual images in each band and 
for each chip. The ~200 stars used to model the point- 
spread functions (PSFs) were carefully selected to cover 
the whole field of view and sample the spatial variations 
of the PSF. The input list of stars for ALLFRAME was 
created with DAOMASTER from the ALLSTAR pho- 
tometry files of the individual images. This list contains 
all the stars that were detected on at least one image. 
The stars with good photometry were selected among 
the detected objects using ALLFRAME's fitting param- 
eters fj, and SHARP. Only those --30000 objects 
with very good photometry, i.e., with a < 0.15, ^ 1-1 
and -0.3 < SHARP < 0.3, were kept. 

Our photometry of chips 1, 3 and 4 of the WFC* was 
calibrated to standard magnitudes using OGLE's pho- 
tometry of the same field, kindly provided by Dr. A. 
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Fig. 1. — Color-magnitude diagram of IC1613. All the stars 
described in section [2] are plotted. 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

See the WFC User Notes at 

|http : // www .ing.iac.es /Astronomy /instruments / wfc / wfc-notes_apr"98! 
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Fig. 2. — Distribution of resolved stars in IC1613. North is up 
and East to the left. 



Udalsky (private communication). Several hundreds of 
stars were used for each chip, giving dispersions of the 
fits at the centers of mass of the point distributions of 
about 0.001. Calibration of chip 2, which was outside 
ogle's field, was more laborious. In V, one standard 
star field observed during the IC1613 INT run was used 
to determine the transformations between chip 2 and 
chips 1 and 4. The dispersion of these transformations is 
'--^0.005. Then the transformation for the latter chips 
based on OGLE's photometry was used, bringing the 
chip 2 V photometry into the standard system. In /, 
chip 2 was calibrated differentially with respect to chip 
4 using overlapping images obtained on the IAC80 tele- 
scope at Izana, Tenerife, Spain, during a photometrical 
observing run. The dispersion at the center of mass is 
of the order of 0.01. He nce, the total error in our pho- 
tometry is that given bv lUdalski eiTall (|200lh . i.e., up to 
0.02 for both V and / bands. Figure [T] shows our final 
{V-I, I) CMD^. The spatial distribution of these stars is 
presented in Fig. O 

The errors given by ALLFRAME are the residuals 
of the PSF-fitting procedure, so they should be con- 
sidered internal errors. Signal-to-noise limitations, stel- 
lar crowding, blending and starloss, which we can refer 
to as observational effects, are important error sources 
and significantly modify the CMP shape an d stellar den- 
sity distribution (jAparicio fc GallaT^ll995f ). To estimate 
the observational effects and the completeness of our 
hotometr y, we re s orted to artificial stars tests. See 



Aparicio fc GallartI ()1995f ) for a detailed description of 
the procedure and the effects of crowding. Basically, a 
large number of artificial stars covering the same range in 
color and magnitude as the observed stars is added to the 
images using the corresponding PSF. These were placed 
on the images following a triangular grid in order to avoid 
crowding between artificial stars themselves and to op- 
timally sample the chip fields. The photometry is then 

^ The photometric data are available from the first author upon 

htnr^ uest. 
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Fig. 3.— Left: Synthetic CMD produced by lAC-STAR showing 
the 153260 stars injected in chip 4 (see text for details). Right: 
CMD of the ^40000 artificial stars recovered both in V and I and 
considered to have good photometry as described in section 2. The 
hne shows the 25% completeness hmit as determined from the ar- 
tificial star tests. 

repeated in the exact same way as was done originally. 
A comparison of the input and output artificial-star list 
gives information about the completeness as a function 
of magnitude and galactocentric distance. 

As the artificial star sam ple, we took a synth etic CMD 
produced by lAC-STAR^ (lAparicio fc Gallartl[20Cil us - 
ing the stellar evolution lilarary of [Bertelh et al 
and t he bolometric corrections of CastcUi & Kurucz 
((200l . The SFR was chosen constant between 13 Gyr 
ago and now, while the metallicity range increases from 
0.0008 to 0.002 at i = 0, to 0.0008 to 0.006 at the present 
time. These ranges were chosen wide enough so that they 
include the actual metallicity of IC1613. In total, 76630 
artificial stars were added in 5 runs in each external chip, 
and twice as many in the central chip (#4) where crowd- 
ing is more important. Figure[3]presents the injected and 
recovered CMDs for chip 4. The broken line indicates the 
25% completeness limit, averaged over the galactocentric 
radius range, obtained from the ratio of recovered to in- 
jected stars as a function of magnitude. Figure H] shows 
the completeness and error in recovered magnitude as a 
function of input magnitude for an inner and an outer 
field. It is important noticing that, even though the re- 
covered artificial stars were filtered using the same values 
of CT, and SHARP as the real stars, external errors 
can be as large as ~1 magnitude at the faint limit. 

3. THE COLOR-MAGNITUDE DIAGRAM 

The most evident feature of our CMD, Fig. [1] is the 
red giant branch ( RGB), composed of low mass stars 
(M < 1.8 - 2.OM0; IChiosi et al]|1992D older than about 
1 Gyr and burning hydrogen in a thin layer around an in- 
ert helium core. The age-metallicity degeneracy and the 
presence of the fainter extension of the asymptotic giant 
branch (AGB) in the same region of the CMD makes it 
difficult to get detailed information about the stars pop- 
ulating it. However, it is possible to place rough limits in 
time for given metallicities using theoretical isochrones. 
A more detailed determination of the chemical enrich- 
ment law (CEL) using this method, as well as its limita- 

^ The code, available for free use, can be executed at the lAC- 
STAR website |http: / /iac-star.iac.es/| 
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Fig. 4. — Completeness (left) and errors in recovered magnitudes 
(right) for a central field (chip 4, top panel), and an outer field 
(chip 1, bottom panel). The error bars in the right panels show the 
dispersion per magnitude bin. 



tions, is presented in section [5. II 

The presence of a well populated main sequence (MS) 
blueward oiV — I = and up to / ~ 18.5 indicates very 
recent star formation (<10 Myr). 

The stars with < {V - I) < 0.6 and / < 22 are most 
likely blue loop (BL) stars, highlighting the blue edge of 
the core He-burning loop, while the red edge (i.e., the 
red supergiant branch, RSG), although contaminated by 
foreground stars, is well defined from (V — I) ^ 1.0 to 
{V — I) ^ 1.8 and up to about / ^ 16.5. These are 
young and intermediate- to high-mass stars, so they are 
among the most metal-rich stars in IC1613. However, 
this is a poorly understood phase in stellar evolution 
and the theoretical models still contain large uncertain- 
ties owing to the importance of processes such as mass- 
loss, overshooting and rotatio n in very massive stars (see 
e.g.. iMaeder fc Mevnet|[200ll) . The BL and RSG actu- 
ally r each the red clump down to / = 23.76 (jCole et al.l 
Il999t ). but the dispersion at the faint end of our CMD 
makes it impossible to distinguish them from the RGB 
stars below / = 22. 

The fi nal extension of the AGB, or red-tail (jCallart et al.l 
Il994r ). extends horizontally redward from the RGB tip at 
/ ^ 20. AGBs are shell H- and He-burning stars, of low 
and intermediate mass and age over about 0.1 Gyr. Their 
relatively large number indicates a possibly important 
intermediate-age population with relatively high metal- 
licity, which would be compatible with the enhanced SFR 
between 3 to 6 Gyrs ago found bv lS03l 

Because of the wide field of view, the diagram also con- 
tains a relatively large number of foreground stars. Most 
stars with 0.5 < V — I < 1.2 and / < 20 are probably 
Galactic dwarfs since their number is the same for the 
middle and bottom panels in Fig. [5] after correcting for 
the difference in area. More generally, the foreground 
contamination is rather important for V — I > 0.6. 

Although old stars (>10 Gyr) are very likely present 
in the RGB, our CMD is not deep enough to detail the 
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Fig. 5. — Color-magnitude diagram of IC1613 for different 
galactocentric distances. Tfie radius ranges are (a) r < 10', 
(b)lO' <r < 16'.5, (c) r> 16.'5. 



old, low metallicity population as it would if fainter MS 
stars were observed. The existence of a bona fide old 
population in IC 1613 is known thanks to the presence o f 
RR Lyrae stars (|Saha et all Il992l : iDolphin et all l2001h 
and CMDs from Hubble Space Telescope observations 
showing core h elium burning, horizontal branch stars 
(|Cole et al. 1999: Dolphin et al. 2001^ , S05) and the old- 
est MS turnoffs (lGallart|l2QQ^ 
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4. MORPHOLOGY, SPATIAL EXTENT & 
DISTRIBUTION OF STELLAR POPULATIONS 



Fig. 6. — Number density of stars as a function of galactocentric 
distance before (triangles) and after (diamonds) background sub- 
traction. The right-hand vertical axis gives a rough estimate of the 
surface magnitude, calibrated as described in text. The horizontal 
line is the weighted mean density of ellipses 29 to 41, which has 
been adopted as the background level. The exponential fit has a 
scale length of 2.'9 it 0.1. The vertical dotted line shows where the 
foreground stars start to dominate, as obtained from the CMDs 
(see section I3J. 



To characterize the morphology of IC1613, we plotted 
the RGB star distribution obtained from our photome- 
try, and convolved the resulting map with a Gaussian 
of o- = 50" using IRAF's GAUSS from the IMFILTER 
package. The result of this process is a smooth map of 
the stellar density highlighting the morphology of the 
galaxy. Fitting ellipses to the isodensity contours was 
done with IRAF's ISOPHOTE routine. 

From the best fitted ellipses, where the crowding is low 
but the star number sufficient, we find a position angle 
of 80° and an eccentricity e = 1 — b /a = .15, in good 
agreement with the values given bv lAblii {PA = 81°, 
e = 0.18; Il972i) . Following the shape of the isopleths, 
we divided the galaxy into concentric ellipses at small 
radii (<10.'5) and circles at larger ones that we used for 
radial star counts and stellar populations gradient analy- 
sis. Their semi-major axis increases in steps of 100 pixels, 
corresponding to 33". 

Figure [6] shows the radial profile of the galaxy con- 
structed from the number of stars in each ellipse after 
correcting for completeness. The correction was obtained 
from the ratio of the number of recovered to injected 
stars in the artificial star test for each annular region. 
The area of the ellipses has been calculated via Monte- 
Carlo sampling, carefully taking into account the gaps 
between the chips as well as the regions around saturated 
stars when calculating the effective surface. An approx- 
imate surface brightness scale, shown on the right-hand 
side of Fig. [SI was calculated from the stellar density in 
each ellipse. The calibration was calculated by compar- 
ing the star number with the total, sky subtracted flux 
for each ellipse then averaged over the radius range for 
consistency. An exponential least square fit to this curve 
between 2' and 15' from the center gives a scale length 
of 2.'9 ±0.1 (620 ± 20 pc) and central surface brightness 
lj,v,o = 22.7 ± 0.6. The former value is sli ghtly smaller 
than t he 760 ± 50 pc (-3f5 ± 0.2) given bv lHodge etlHI 
()1991h . However, a similar larger length is obtained when 
fitting only the inner ^7', as was done by the authors 
because of the limited depth of their photometry, and 
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Fig. 7. — Spatial distribution of young (MS+BL+RSG; top) 
and older (AGB+RGB; bottom) stars. The red ellipse shows the 
extent of the MS stars (r=10'). The red dashed circle indicates 
where the foreground stars start to dominate (r=16.'5). The neutral 
hydrogen emis s ion co ntours (0.2, 6.2 and 14.2 Jy km s~^) from 
I Hoffman et al.l l|1996l ) have been overplotted in blue (see text for 
details). 



omitting the crowding correction. The profile seems to 
get steeper at r 12', but this change of slope could 
be an artifact of the background subtraction and small 
number statistics. 

The large field of view covered by the WFC per- 
mits to study the spatial gradients of the stellar pop- 
ulation across the galaxy. Ideally, that would give 
us hints on its formation and evolution. However, in 
the case of shallow CMDs, the spatial variations of 
the morphology of the CMD only reflect accurately 
differences in the SFH over the last several hundred 
million years. Such variations have been observed in 
all the known dwarf irr egular galaxies through CMD 
morphology (e.g., WLM: iMiWiti fc Ziilstral[l996l : NGC 
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Fig. 8. — Stellar densities of the different populations of resolved 
stars after correction for completeness and background contami- 
nation. The stars selected for each population are shown in the 
inset. The RGB profile was divided by four to fit on the graph. 
The errorbars have been omitted for clarity. Note the difference 
in the scale length of the older (RGB, AGB) and younger (MS, 
BL-(-RSG) populations (see text for i 



6822 : iBattinelli et~al] 120061: Leo A: IVansevicius et all 
12004 Phoenix: iMartmez-Delgado et al.l 119990 or distrT 
bution of the varia, ble star p opulations (e.g., in P hoenix: 
iGallart et ai]|200l in Leo L iBaldacci etal]l200l . 



Figure [H] shows that it is also the case for IC1613 and 
confirms the diffe rence in relative number of young and 
old stars found bv lS03l between their outer H ST field and 
a central field studied earlier by Cole et all (jT999). 



We 

divided the galaxy into three regions following the mor- 
phology of the CMD of each individual ellipse defined 
above: the inner part of the galaxy where stars younger 
than about 500 Myr are present (r < lOfl, Fig. [5^), the 
region at intermediate distance with no young star but 
still a well defined RGB (lOfl < r < 16.'5, Fig. [Sb), 
and the outermost part of our observed field, dominated 
by foreground stars [r > 16.'5, Fig. [SJ;). The corre- 
sponding spatial limits, as well as the HI contours from 
[Hoffman ct al. (1996), are displayed over the stellar dis- 
tribution of the young and old populations in Fig. [71 
Note that the HI contours were shifted by — 0.'6 (-1-100 
pixels in y) and 2'.9 (—500 pixels in x) in right ascen- 
sion and declination, respectively. Th is c orresponds to 
the of fset between the as t romet ry of [Lake & Skillma^ 
(1989) and [Hoffman et al.l ([1996f ). and adopting the for- 
m er for the H I contou rs correctly places the H II regions 
of [Hodge et al.[ ([1990| ) on top of the star forming regions, 
and fits our star distribution better. Although the RGB 
of Fig. [St is not clearly defined, a substantial fraction of 
the displayed stars prob ably belong to IC16 13. This cor- 
roborates the results of [Albert et al.[ ([20001 ) . who found 
carbon stars extending up to 15' from the center of the 
galaxy, and shows that IC1613 is actually more extended 
than previously thought. 

To give a quantitative measure of the gradient in the 
age composition of the stars in IC1613, the stellar sur- 
face density for different populations of resolved stars is 
presented in Fig. [SI The age gradient is clearly visible: 
while the density of the older stars follows the expected 
exponential decrease from the central region, that of the 
young population peaks at a radius of ~3' and vanishes 
rapidly as the radius increases. This results in the scale 



6 



Bernard et al. 




0.25 r 



0.20 r T 

12 3 4 

T [Gyr] 

Fig. 9.— Top: Width of the RGB of IC1613 at Mj = -3.5 as 
a function of galactocentric distance, before {thin line) and after 
{thick red line) correction for the effects of crowding. The dashed 
hne indicates the mean value of the corrected width beyond ~5'. 
Middle: Width of the injected {long-dashed line) and recovered 
{full thin line) RGB of the artificial stars tests. The thick red line 
is a fit to the recovered width. Bottom: Width of an artificial 
RGB as a function of the age of the youngest stars in the CMD. 
The small change between t = and 1 Gyr is probably due to the 
presence of stars from the RSG sequence. In all cases the errorbars 
correspond to the dispersion of twenty solutions. The vertical scale 
is the same on all three panels. 

length of the young population being much smaller than 
that of the older stars: a fit to the profiles between 3f5 
and 12.'1 from the center of IC1613 gives scale lengths of 
1.'19 ± 0.04, 2f0 ± 0.1, 3.'2 ± 0.2 and 3f8 ± 0.1 for the MS, 
BL+RSG, AGE and RGB, respectively. The off-centered 
peak in the distribution of the young population, visible 
at r ~ 3'-3.'5 in Fig. [8l is due to the fact that the star 
forming regions are distributed in a some what circular 
pattern at this di stance from the center (Hodge et al.' 
[1990), where Silic h et al.l ([2006) observed a higher H I 
column-density. 

Because of the large dependence of the RGB color on 
metallicity, the position of the RGB in a CMD can be 
used to estimate the metallici ty of the stell ar system 
using empirical relations (e.g., iDa Costa fc A rmandrofj 
|1990| ). For the same reason, the width of the RGB is 
often considered to be an indication of the metallicity 
dispersion. On the CMDs presented in Fig. [5l it appears 
that the RGB in the central region is wider than that 
at larger radii and suggests the existence of a metallic- 
ity gradient. However, a RGB width gradient may also 
reflect an observational-effects gradient. To check this 
possibility, the following procedure was used, including 
measurements of real- and artificial-star RGB width and 
error estimates. We measured the width of the RGB at 
Ml = —3.5 as a function of galactocentric distance for 
the real stars as follows: for each radius interval, cho- 
sen so that the RGB contains about 800 stars, a sub- 



sample of 450 RGB stars was selected at random. The 
resulting RGB was sliced in intervals of 0.25 magnitude 
between 20.5 < / < 22, and the width of the RGB in 
each magnitude range was obtained from the full-width 
at half-maximum (FWHM) of a Gaussian fit to the color 
function (CF). The width at Mi = —3.5 was then inter- 
polated from a linear fit to the FWHM versus magnitude 
plot. We repeated this operation twenty times using a 
different subsample of RGB stars each time, and used 
the dispersion of these solutions as the errorbars. The 
same process was followed for the synthetic CMD used 
in the artificial stars tests, except that the radius bins 
and random subsamples contained 2500 and 1500 stars, 
respectively. The resulting plot is presented in Fig. [5] for 
IC1613 {top) and the artificial stars [middle panel). 

It shows that the RGB recovered in the artificial star 
tests is wider in the center of the galaxy. In addition 
to the expected FWHM widening due to signal-to-noise 
limitations, crowding in the central part further disperses 
the stars on the CMD. Crowding is important and af- 
fects the width of the RGB up to r ~ 10'. To correct 
the observed RGB for the effects of crowding, we calcu- 
lated the broadening parameter as a function of radius 
by subtracting in quadrature the injected width from the 
recovered one. The corrected width was then obtained 
by subtracting in quadrature the broadening parameter 
from the observed RGB width, and is shown as the thick 
red line in the upper panel of Fig. O The FWHM of 
the corrected RGB still presents a significant variation 
across the central '^6'. A higher metallicity dispersion in 
the central region could be responsible for this residual 
width excess. However, the lower panel of Fig. [9l show- 
ing the width of artificial RGBs for which the age of the 
youngest stars is progressively older, indicates that this 
excess is consistent with the presence of stars between 1 
and 3-4 Gyr old in the center. Thus, there is no strong 
evidence of a metallicity gradient across the galaxy but 
its presence cannot be ruled out with the present set of 
data. 

5. RECENT STAR FORMATION HISTORY 
5.1. The Method 

The reconstruction of the recent SFH has been per- 
formed in a way similar to that described in Hidalgo & 
Aparicio (2007, in preparation) or Gallart et al. (1999). 
The main difference here resides in the fact that our 
CMD is not sufficiently deep to contain discriminating in- 
formation about intermediate- age and old stars. We lim- 
ited our study to the recent SFH, i.e., the last ~300 Myr, 
based on the MS, BL and RSG populations. It thus lim- 
ited the determination of the SFH to the inner r < 10' 
since these populations are absent beyond this radius. 

In short, the SFH is derived from the comparison of 
the star distribution of the observed CMD with that of 
a synthetic CMD. The synthetic CMD to be used in the 
comparison needs to cover a range of age and metallic- 
ity at least as large as the one that can be expected 
in such a d warf irregular galaxy. It was generated by 
lAC-STAR dAearicio & Gallart' '2004") using the stellar 
evolution library of Bertelli et al. (1994) and the bolo- 
mctric corrections from iCastelli fc Kuruc^ (|2003l ) , with 
the following parameters: the SFR was chosen constant 
between 13 Gyr ago and now, the initial mass function 
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Fig. 10. — CEL of IC1613. The diamonds indicate, for a given 
age, the minimum and maximum values found by fittin g theoretical 
isochrones from the library of Padua IjGirardi et al.ll200a ) to the 
upper RGB. For reference, the CELs obtained via the 'Tolstoy' 
and 'Dolphin methods' from |S0 3| are shown as red circles and blue 
triangles, respectively. The grayed region is the CEL used for the 
comparison CMD in the SFH analysis. The stripped region in the 
last 300 Myr shows the metallicity range that gave the lowest 
when determining the SFH. The right-hand scale was converted 
from the abundance by mass assuming Zq = 0.0198. 

(IMF) was that of lKroupa et al.l ()1993[ ). and the fraction 
of binaries was set to zero. To fix the input CEL we 
fitted isochrones of different metallicities to the upper 
RGB of the o uter field whe re crowding is not dominant. 
We used the iGirardi et al.l (2002) isochrones for differ- 
ent ages (1.585, 2.239, 5.012, 10.00 and 12.59 Gyr old) 
and metallicities between Z = 0.0005 and Z = 0.006 
in steps of 0.0005. The metallicities that are not avail- 
able in the original library were interpolated using lAC- 
STAR. A metallicity was considered valid at a given age 
if the corresponding isochrone was inside the FWHM of 
the RGB between 20.4 < / < 21. Figure [10] (/^/ted dia- 
monds) shows the values obtained through this method. 
Although the method is rather simplistic, our resulting 
CEL is in very good agr e emen t with those derived by SQ3- 
However, ICallart et al.l ()2005f ) showed that the theoreti- 
cal isochrones are generally too vertical with respect to 
the empirical ones, leading to a slight overestimation of 
the metallicity. The grayed region represents the metal- 
licity range employed to create the comparison CMD. 
It is mainly used for producing the RGB, from which 
we obtain the mean SFR between 1500 and 13000 Myr, 
and serves a normalization purpose for the young SFH. 
For the stars younger than 300 Myr, we tried to put 
constraints on their metallicity by further restricting the 
metallicity range in the comparison CMD to intervals of 
0.002 in abundance by mass, and ran the algorithm three 
times with the following ranges: Z = 0.0015-0.0035, 
Z = 0.0025-0.0045 and Z = 0.0035-0.0055. 

In order to simulate the observational effects and al- 
low a more realistic comparison with the real CMD, we 
applied the dispersion in color and magnitude of the syn- 
thetic stars recovered from the crowding tests to the 
comparison C MD, following the same procedure as in 
ICallart et al] (1999). 

In the observed and comparison CMDs, the MS, BL 
and RSG sequence are then divided into 'boxes'. Because 
the box selection could influence the resulting SFH, we 



used three different types of box to rule out this possi- 
bility: two regular grids with large (0.3 and 0.5 in color 
and magnitude, respectively) and small (0.2 and 0.3) box 
size, and an "a la carte" parametrization. In the lat- 
ter, the size and shape of the boxes are chosen taking 
into account the knowledge and model limitations of the 
evolutionary phases, in particular the slope of the RSG 
sequence and the position of the BL. The regular grids 
were also shifted in color and magnitude to check for 
consistency and ensure the significance of the result. 

Additionally, the synthetic stars are divided into par- 
tial models, each with a small range in age. The tempo- 
ral resolution is limited by the quantity of information 
contained in the CMD, which depends on the quality of 
the data. It was chosen by comparing the capacity of 
the algorithm to recover the known SFH of a synthetic 
CMD, after the observational effects had been simulated 
as described above, using different time ranges. Above a 
certain time resolution — which depends on age — the re- 
covered SFH was mainly made of short, violent bursts 
of star formation separated by peri ods of inactivity, re- 
gardl ess of the input SFH (see also lAparicio fc Hidakol 
I200I . 

The reconstructed SFH is a linear combination of the 
different partial models. The best solution is obtained 
through x^-minimization by a genetic code. A thorough 
description of the algorith m and method is presented in 
lAparicio fc Hidalgol (|2007D . In total, we used 24 models 
with different CMD parametrization and time resolution. 
The consistency of the method was checked by solving 
the SFH of synthetic CMDs and comparing the solutions 
with the input SFHs. Some of the regular grid models 
did not give a satisfactory solution, but the discrepan- 
cies could generally be traced to small differences in the 
location of evolutionary phases in the CMD, while the 
solutions obtained with the d la carte parametrization 
were found to be more stable and are the ones we will 
present here. 

5.2. The Results 

Of the three metallicity ranges used for the stars 
younger than 300 Myr, the second one, i.e., Z = 0.0025- 
0.0045, gave the best solutions, and is shown as a 
stripped region in Fig. ITOl This is co nsistent with the 
H II region metallicity of lTalentI (|1980l ). Our best solu- 
tions for the SFHs at different galactocentric distances 
are presented in Fig. 1111 The radius ranges were cho- 
sen so that the CMDs contain the same number of stars. 
The error bars correspond to the dis persion of twenty so- 
lution s for which = xlest + ^ f see lAparicio fc Hidalgol 
[2OOI . 

The overall picture is a relatively constant SFR at all 
radii, decreasing with increasing radius, while the mean 
age of the stars increases with radius. The very central 
region, panel (a), shows a fairly constant SFR for the 
last 300 Myr of (1.6 ± 0.8) x lO'^ M© yr"! kec^!^ 
excellent agreement with t he va lue found by ICole et al.l 
(1.6x10-3 Mq yr-i kpc-^: ll999f ) using the V-hand lumi- 
nosity function of their WFPC2 central field and assum- 
ing a Salpeter IMF. The sharp drop in the last ~25 Myr 
is in agreement with the lack o f very bright stars at the 
center of IC1613 first noted bv lHodge et all 11991). 

Between 2.'5 and ~6' [(b)-(d)] the SFR increased by a 
factor ~2-3 in the last 100-150 Myr. This corresponds 
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Fig. 11. — Best SFR(t) obtained for each galactocentric distance. The radial ranges were chosen so that the corresponding CMDs contain 
the same number of stars. The error bars correspond to the dispersion of twenty solutions for which Xa = + 1) where is the residual 
of the best solution, indicated in each panel. The gap in the solutions between 300 and ^1500 Myr i s due to th e lack of information f rom 
either the young populations or the RGB in this age range. The solutions obtained by ICole et al.l | |T999D and ISkillman etH] 112003) for 
their respective WFPC2 field are plotted in panels (a) and (e). 



to the peak in the ra dial distributi o n of y oung stars and 
HII regions found by [Hodge et al.l (|199a) . 

The field studied by lS03l is a small fraction of the region 
for which the SFH is represented in p anel (e). The SFH 
calculated via the 'Cole method'^ in IS03l . shown as the 
green long-dashed lines, is very similar to the one we 
obtained here for the whole elliptical annulus. At larger 

''' Three methods were used in IS03I to calculate the SFH: the 
'Dolphin method' gives the SFR on a relative scale only. For the 
'Tolstoy method', the vertical scale of their Fig. 7 gives a mean 
SFR ~ 3 X 10-3 M0 yr-l kpc-2, a factor ~15-20 higher than 
that of the 'Cole method' and as high as the mean value of the 
SFR over the whole galaxy IIMatedll998l '). 



radii, the number of stars formed more recently than 
about 300 Myr decreases to a negligible value when the 
radius reaches ^10'. 

6. CONCLUSIONS 

We have presented an analysis of the stellar content, 
morphology, and recent star formation history of the Lo- 
cal Group dirr galaxy IC1613 based on wide-field {V-I,I) 
photometry of resolved stars. 

The distribution of resolved stars can be fitted with 
ellipses of position angle 80° and eccentricity 0.15. The 
exponential fit of the resulting radial profile has a scale 
length of 2.'9 ±0.1. The relatively large number of RGB 
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stars still present in the outer part of our observed field 
(r > 16f5) indicates that the galaxy is actually more 
extended than previously estimated. 

The well-populated young evolutionary phases of the 
CMD of the central region of IC1613 indicate very re- 
cent star formation (<10 Myr). The changing CMDs as 
a function of galactoccntric distance show a strong gradi- 
ent in the age of the younger stellar population, with the 
young stars lying preferentially in the central part, while 
the older ones are distributed more uniformly. No evi- 
dence of star formation more recent than about 300 Myr 
was detected beyond r > 10'. 

Analysis of the width of the RGB as a function of ra- 
dius is consistent with no metallicity gradient. The com- 
bination of crowding effect and the presence of younger 
stars in the RGB is responsible for the widening toward 
the center of the galaxy. 

In the region where the recent SFH could be studied 
(r < 10'), the results indicate a decreasing SFR(t) from 



the center outward as expected from the distribution of 
neutral gas in IC1613, with the exception of the annular 
region where the star forming regions are clustered (r ~ 
4') and therefore the SFR is a factor ~2-3 higher. 
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